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Abstract. We discuss the science motivations and prospects for a joint analysis of 
gravitational-wave (GW) and low-energy neutrino data to search for prompt signals 
from nearby supernovae (SNe). Both gravitational- wave and low-energy neutrinos are 
expected to be produced in the innermost region of a core-collapse supernova, and 
a search for coincident signals would probe the processes which power a supernova 
explosion. It is estimated that the current generation of neutrino and gravitational- 
wave detectors would be sensitive to Galactic core-collapse supernovae, and would also 
be able to detect electromagnetically dark SNe. A joint GW-neutrino search would 
enable improvements to searches by way of lower detection thresholds, larger distance 
range, better live-time coverage by a network of GW and neutrino detectors, and 
increased significance of candidate detections. A close collaboration between the GW 
and neutrino communities for such a search will thus go far toward realizing a much 
sought-after astrophysics goal of detecting the next nearby supernova. 
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1. Motivation 

The predicted rate of core-collapse supernovae (SNe) in our Galaxy is ~ 2 per century; 
the rate is about twice the Galactic rate out to the Andromeda galaxy (~ 1 Mpc), 
and about one per year out to the Virgo cluster (~ 10 Mpc) pp. These numbers are 
accompanied by large uncertainties inherent in converting galaxy properties to supernova 
rates. However, [1] also suggests that the increased number of nearby core-collapse 
supernovae discovered within the past few years strongly indicates that the predicted 
rates might be significantly underestimated by a factor of ~ 3 in the 3 — 5 Mpc range. A 
direct upper limit on the Galactic SNe rate, based on non-observations of antineutrino 
events in the past 25 years, is given in [2]. 

Contemporary neutrino detectors and kilometer-scale gravitational-wave (GW) 
detectors are currently poised to detect the next Galactic core-collapse supernova. Both 
neutrino and gravitational-wave signals are expected to be generated in the innermost 
region of a dying star, and both signals are expected to be emitted within a short 
time interval of each other, i.e. within a few milliseconds [31 H]. While both neutrino 
and GW detectors are preparing to independently detect a Galactic supernova, there 
are science benefits to systematically searching for a supernova signature using a joint 
analysis of neutrino and GW data which is guided by the expected proximity of the 
neutrino and GW signals. These science benefits include lower detection threshold 
requirements, better live-time coverage, increased significance of candidate detections, 
extended distance reach to the local volume of galaxies, and increased sensitivity to 
core-collapse events which have only a weak or non-existent electromagnetic signature. 

The detection of a burst of low-energy neutrinos from SN1987A, at a distance of 
about 50 kpc in the Large Magellanic Cloud, by the Kamiokande II and Irvine-Michigan- 
Brookhaven (1MB) experiments, and by scintillation neutrino detectors [51 El El |9] 
demonstrated the capability of neutrino detectors to detect SN events, and paved the 
way for a description and validation of the standard model of neutrino emission from 
a core-collapse supernova. In a core-collapse supernova, ~ 99% of the neutron star's 
gravitational binding energy (~ 3 x 10 53 ergs) is released in the form of neutrinos (and 
anti-neutrinos) of all flavors. These neutrinos have energies in the few tens of MeV 
and are emitted over a time scale of a few tens of seconds. The neutrino light curve is 
expected to show structure, with an increase in luminosity during the first ~ 0.5 second 
due to accretion of matter onto the proto-neutron star [TOl \TT\ [12] . A neutronization 
burst, which is a peak in the v e luminosity produced as the shock from the core bounce 
propagates through the star's outer core, is expected to last a few milliseconds after the 
core bounce. However, the energy emitted in the neutronization burst is only < 1% of 
the total energy emitted, and such a feature might not be easily recognized in a measured 
neutrino light curve because of its short duration, and because the cross section of v e e 
scattering is lower than that of the dominant inverse beta decay (u e p — > e + n) reaction 
in a neutrino detector's medium [10]. The events detected by Kamiokande II and 1MB 
from SN1987A are reproduced in figure [TJ 
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Figure 1. The energy of the events detected by Kamiokande II (closed circles) and 
1MB (open circles), produced by the neutrino burst from SN1987A, as a function of 
time, in seconds [7]. 



The groundbreaking discovery of a neutrino burst from SN1987A was guided by an 
optical sighting of the supernova [3 [7]. The optical sighting of such a close astrophysical 
event motivated the analysis of archived neutrino data and guided the time scale in 
which to perform the search. It is plausible, however, that a fraction of core-collapse 
events are accompanied only by a faint electromagnetic display. This might be due 
to extinction brought about by an extremely dusty environment or the intervening 
interstellar medium, or an inherently weak accompanying electromagnetic emission with 
a fast decay time p2]. For example, the supernova of Cassiopeia A, one of the youngest 
known Galactic supernova remnants and which is at a relatively close distance of 3.4 kpc, 
has no historical record of widespread sighting [14j [15] during the epoch when the 
explosive fireworks would have reached Earth around 340 years ago, and it is plausible 
that this supernova fell in this category [TU [161 02] • Obscuration by dust could also 
explain the non-sighting of the supernova of an even younger, < 150 year-old remnant, 
Gl.9+0.3 |18j . This would be analogous to how the dearth of known Galactic supernova 
remnants could possibly be attributed to low surface brightness, leaving faint supernova 
remnants unresolved from the Galactic background emission [19] . 

It has also been argued that observations suggest a deficit of optically observed 
high-mass (> 25M Q ) core-collapse SN progenitors [15] . and that optical searches have 
provided little information on the possibility that massive stars end their lives by forming 
black holes without the dramatic electromagnetic signature of an explosion. Reference 
[T3] estimates that an upper bound on the rate of so-called "failed" supernovae is roughly 
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equal to the rate of successful SNe, and that the lower bound on black hole formation 
rate is ~ 25% that of normal SNe. 

A joint GW-neutrino search for nearby core-collapse supernovae could potentially 
provide insight to this scenario where a considerable fraction of stars end their lives with 
little or no electromagnetic display 

Indeed, the natural progression of the expected supernova signature — from the 
prompt GW and neutrino signals to the optical signal which is expected to rise 
many minutes to hours later — and the current state of neutrino and gravitational-wave 
detectors and their respective detection algorithms, make it likely that the detection 
of the next nearby supernova will proceed in a direction opposite that of SN1987A, 

1. e. that prompt neutrino and GW detections would trigger optical telescopes to search 
for an optical counterpart. The infrastracture of the Supernova Early Warning System 
(SNEWS), for example, is designed to alert observatories in the event of a detection 
of neutrinos from a nearby supernova [20l |2T] . Simulations also indicate that, for a 
Galactic supernova, the neutrino detectors Super-Kamiokande and IceCube would be 
able to reconstruct the SN bounce time to within a few milliseconds j22j [23]. On the 
gravitational- wave side, there is also an alert system, called LOOC UP [21], which is 
being developed to send plausible future candidate GW triggers to optical observatories 
for confirmation of a corresponding astrophysical source. 

Astrophysical events — such as gamma-ray bursts (GRBs) and flares from soft 
gamma repeaters (SGRs) — detected by other observatories have been extensively 
utilized as external triggers in the analysis of LIGO- Virgo data to search for GW 
counterparts to these events [25j [2HJ [23 [23 123 SO]. The strategy of using external 
triggers with precise event timing and position information to look for GW signals is 
motivated mainly by the decrease in both background rate and the effective number of 
experimental trials that shorter analysis time windows make possible. In the case of 
GRBs and SGRs, the analysis windows range from a few seconds to a few minutes. On 
the other hand, the use of an optical signal from a supernova as an external trigger does 
not provide the same tight constraints on the time and duration of the analysis window. 
Studies indicate that, at best, the time of a supernova explosion can be determined from 
an optical light curve to within a few hours, but only if the first measurement of the 
optical flux is made within a day of the explosion |31j . 

In contrast, the gravitational- wave and neutrino signals are expected to be detected 
within a tight window, ranging from a few milliseconds to a few seconds, depending on 
the dominant GW emission process [31 E]. 

2. Science benefits of a joint GW-neutrino search 

Several of the world's neutrino detectors have performed searches for core-collapse 
supernovae and have evaluated their respective detection efficiencies as a function of 
distance. The Super-Kamiokande (Super-K) water Cherenkov detector in Japan [32J, 
the scintillation detectors Large Volume Detector (LVD) |33j and Borexino [31] in Italy, 
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and KamLAND in Japan [35], IceCube at the South Pole [36], MiniBooNE in the USA 
[37] . and others [38] are expected to robustly detect a neutrino burst from a Galactic 
supernova in the baseline model scenario. Super-K, for example, would detect ~ 8000 
events for a core-collapse SN at the center of the Milky Way, ~ 8.5 kpc away [32] . The 
Baksan scintillation detector had previously also performed a systematic search [39]. 

Analogously, all-sky searches for GW bursts have been performed using data from 
the Laser Interferometer Gravitational- Wave Observatory (LIGO), the most recent 
of which made use of data from the fourth (S4) and fifth (S5) LIGO science runs 
[4"0~1 |4"T] . While no gravitational waves have been directly detected from an astrophysical 
source, the current generation of LIGO-Virgo interferometric GW detectors have made 
tremendous progress in improving their sensitivities [121 H31 HH 051 H6], and are 
expected to be sensitive to several models of GW emission from a Galactic core-collapse 
supernova [3111]. Significant improvements in sensitivy are expected to continue with the 
anticipated advent of the next generation of GW interferometers, Advancd LIGO and 
Virgo. However, while the astrophysical motivation for expecting gravitational waves 
to accompany core-collapse supernovae is strong, the expected rate, gravitational-wave 
strength and waveform morphology are uncertain [3111]. As a benchmark, the expected 
energy going into gravitational waves may range from 10~ 10 to 10~ 4 M c 2 (or 2 x 10 44 
to 2 x 10 50 ergs), and thus only a small fraction of the energy liberated in neutrinos. 

Estimating the sensitivity improvement of a gravitational- wave search due to a 
tighter search window and lower background rate requires assumptions on the spectrum 
of the background events. In the all-sky GW burst search using the first-year data of 
LIGO's fifth science run, a false alarm rate of 1 in 100 years in the 64-200 Hz frequency 
band corresponds to a three-detector network signal-to-noise ratio (SNR) threshold of 
~ 8.5. For frequencies above 200 Hz, the corresponding SNR threshold is lower, but 
the interferometers' strain sensitivity is lower at this frequency band (cf. figure 5 and 
figure 2 of [41], and also Appendix E of [4T]). Requiring a coincidence of GW events 
within O(ls) of a neutrino signal tuned at a rate of 1 per day would allow GW detectors 
to operate at a false alarm rate of 3 x 10~ 5 Hz, which in turn corresponds to a SNR 
threshold of ~ 3.5 [H], or an improvement of a factor of ~ 2 in sensitivity. The 
distance reach of gravitational-wave detectors scales linearly with the inverse of SNR. 
Such potential improvement in gravitational-wave sensitivity, in a joint GW-neutrino 
search, will increase the science reach of the GW instruments relative to what they can 
achieve alone. 

To quantify the science reach, it is important to appreciate that the guidance from 
source phenomenology is subject to significant uncertainties [4J. All estimates of GW 
bursts associated with supernovae rely on models (see [171 HH HH [501 EH E21 [53] for 
recent GW emission estimates). Most such models are not yet 3-dimensional, do not 
incorporate the entire set of possibly relevant physics, and, most importantly, do not 
(in most cases) predict robust supernova explosions as observed in the electromagnetic 
universe. Moreover, none of the current "state-of-the-art" simulations can make reliable 
predictions of the mechanism responsible for the observed velocities of pulsars of 
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Figure 2. Sensitivity of gravitational- wave detectors to core-collapse supernovae as 
limited by the efficiency of past, present and foreseen searches. The diagonal lines 
in this log-log plot reflect the fundamental relation connecting energy going into 
gravitational waves, distance to the source, and root-sum-square strain sensitivity of a 
search, h ISS (cf. equation [TJ . These lines correspond to a fixed search sensitivity at the 
gravitational-wave detectors for narrow-band signals in the most sensitive frequency 
region (~ f50 Hz). All combinations of energy-distance above and to the left of 
these diagonal lines can be probed in a search. The dotted line identified as "S4" 
corresponds to LIGO's fourth science run, and the solid line labelled "S5" to LIGO's 
fifth science run. The dashed line reflects the estimated improvement in sensitivity in 
a joint GW-neutrino search. The dot-dashed line is the expected reach of Advanced 
LIGO and Virgo, with no assumption of joint searches with neutrinos [43j . The 
horizontal lines represent upper bounds on the energy release for four core collapse 
models (as summarized in [3]): A: PNS pulsations (acoustic mechanism, [37]): B: 
rotational instability; C: rotating collapse and bounce; D: convection and SASI. 

up to 1000 km/s (but see [M] who do predict such kicks, albeit with a simplified 
model). It is likely that these velocities were imparted on the neutron stars at birth 
("pulsar birth kicks") which obviously must involve multi-dimensional dynamics and 
gravitational-wave emission presently not accounted for in models. Thus, despite the 
availability of multiple potential explosion mechanisms and their associated multi- 
dimensional dynamics and gravitational-wave signatures, the current picture is unlikely 
to be complete. Moreover, the current, most pessimistic estimates are probably 
overly so in predicting the gravitational-wave yield of core-collapse supernovae. In 
the absence of complete models, observations can and must guide our understanding 
of the astrophysical systems. Joint analysis of neutrino data with initial and enhanced 
LIGO-Virgo observations would significantly enhance the capability to constrain models 
of gravitational-wave emission in core-collapse events. 
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Figure 3. Estimated probability of satisfying a Super-K neutrino burst search 
criterion as a function of distance. Solid curve: standard search parameters 32 . 
Dashed curve: probability if only a single neutrino event is required, in coincidence 
with a GW signal. 



For isotropic emission of gravitational waves, the luminosity distance is related to 
the energy emitted by the source in GW waves, £?gw> and to the root-sum-square strain 
amplitude at the detector, h rss , by [40J, 

£gw^^^ 2 /X s , (1) 

where f is the observed frequency of the waves. For a hypothetical source at a Galactic 
distance of 10 kpc and an assumed signal morphology of a sine-Gaussian waveform [4T] 
with central frequency 153 Hz and quality factor Q = 9 (where Q is a dimensionless 
quantity which is roughly a measure of the width of the waveform in terms of number of 
cycles of the sinusoid), the mass sensitivity during the LIGO S5 run is 1.9 x 10~ 8 M c 2 
|41j . The proposed joint search, with a factor of ~ 2 improvement in sensitivity, would 
decrease by a factor of ~ 4 the minimum energy which is probed in core-collapse 
supernovae, i.e. down to 4 x 10 -9 M Q c 2 for a Galactic supernova with signal content 
in the most sensitive band of the gravitational-wave detectors. Searches at higher 
frequencies would be penalized by the / 2 dependence and by lower strain sensitivity 
of the detectors at these frequencies. At the same time, however, such searches are 
characterized by a lower accidental background rate [41] . and the improvement in 
sensitivity which could be achieved at higher frequencies would be about the same 
as that which was estimated for searches at the instruments' most sensitive frequency, 
i.e. a factor ~ 2 improvement (cf. figure 5 of [H]). These scaling laws, and the 
potential improvement in science reach offered by a joint search, are summarized in 
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figure [21 together with the expected gravitational- wave emission in four sample emision 
mechanisms. 

A gravitational- wave coincidence requirement also has the potential to improve 
the sensitivity of neutrino experiments by relaxing the criteria for detection. For 
example, Super-K's recent "distant" burst search [32J requires two neutrino events (with 
energy threshold 17 MeV) within 20 seconds, which corresponds to approximately 8% 
probability of detecting a supernova in Andromeda. The accidental coincidence rate for 
this criterion is less than one per year; the single event rate at this threshold is about 
1 per day. If one could achieve an acceptable accidental rate by requiring coincidence 
of a single neutrino event with a gravitational- wave signal, then the probability of a 
core-collapse event in Andromeda satisfying the search criterion would be about 35%, 
as shown in figure [3j Distant burst search parameters could be re-optimized with respect 
to current ones; the neutrino event energy threshold could potentially be reduced, further 
improving sensitivity. 

3. Summary 

We have motivated the search for nearby core-collapse supernovae using a joint analysis 
of low-energy neutrino and gravitational- wave data, and we have shown examples of the 
science benefits of such a joint analysis. Turning this idea into a reality in the immediate 
future using contemporary neutrino and gravitational-wave data would make possible a 
richer exploration of the innermost, dynamical processes in a core-collapse supernova. 
A search like this is a necessary complement to the joint high-energy neutrino and 
gravitational- wave search that is currently being planned [55]. Moreover, embarking 
on such a task now would be forward-looking, since this kind of analysis would gain 
importance as the sensitivities of experiments improve. The Advanced LIGO and Virgo 
detectors [56] are expected to start operating in 2014, and are designed to improve on the 
sensitivity of the initial detector configurations by a factor of ~10. At the same time, 
there are a number of large future neutrino experiments planned, employing various 
technologies, including some of megaton scale [571 EH1 EH1 EH EI]. Reference pQ points 
out that such detectors will be able to observe on the order of one supernova neutrino 
event every few years from beyond the Local Group of galaxies (~ 2 Mpc). In such 
a regime, some kind of external (non-neutrino) trigger will be essential to distinguish 
supernova neutrino-induced events from background. 
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